Structural Transitions and Soft Modes in Frustrated DNA Crystals 
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Relying on symmetry considerations appropriate for helical biopolymers such as DNA and fil- 
amentous actin, we argue that crystalline packings of mutually repulsive helical macromolecules 
fall principally into two categories: unfrustrated (hexagonal) and frustrated (rhombohedral). For 
both cases, we construct the Landau- Ginzburg free energy for the 2D columnar-hexagonal to 3D 
crystalline phase transition, including the coupling between molecular displacements along biopoly- 
mer backbone to displacements in the plane of hexagonal order. We focus on the distinct elastic 
properties that emerge upon crystallization of helical arrays due to this coupling. Specifically, we 
demonstrate that frustrated states universally exhibit a highly anisotropic in-plane elastic response, 
characterized by an especially soft compliance to simple-shear deformations and a comparatively 
large resistance to those deformations that carry the array from the low- to high-density crystalline 
states of DNA. 
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Even ignoring its fundamental role as the carrier of ge- 
netic information for many living organisms, DNA is a 
remarkable molecule, which continues to draw significant 
study purely as a component of molecular systems with 
complex order. Although negatively charged in aqueous 
solution, DNA can be condensed by a variety of meth- 
ods P, [4] yielding a remarkably rich spectrum of liquid- 
crystalline order. Roughly speaking, as the concentration 
of condensed DNA is increases a sequence of thermody- 
namic phases can be observed: isotropic — > blue (or pre- 
cholesteric) phases — > cholesteric — > columnar hexago- 
nal — > crystalline 0, 0, [E]. Of course, it is this latter 
crystalline state that is the focus of studies which re- 
solve the double-helical structure of DNA on the atomic 
scale 

The precise nature of the crystalline ordering of DNA 
is quite a complex issue stemming from the microscopic 
frustration introduced by the hexagonal packing of helical 
molecules. While single helices can be uniformly packed, 
the perfect packing of symmetric bi-helical molecules, 
such as F-actin, cannot be achieved in a hexagonal lat- 
tice. Molecular DNA is intermediate to these limiting 
cases although experimental observations of crystalline 
states suggest that interactions between DNA are highly 
frustrated. It is the aim of this Letter to demonstrate 
that as a direct result of this molecular-level packing 
frustration, the 3D crystalline state exhibits an unusual 
response to shear in the plane of hexagonal order, char- 
acterized by a generic softness to simple-shear deforma- 
tions. 

At the high packing densities of crystalline DNA, the 
microscopic interactions between neighboring helices are 
complex, involving a combination of electrostatic, hydro- 
gen bonding, steric and ion-mediated forces. Rather than 
relying on a detailed description of interactions at the 
molecular scale, we appeal to symmetry-based consider- 
ations to explore the consequences of frustration of helical 
molecules. To this end, we construct an order-parameter 
expansion of the free energy near to the 2D columnar- 
hexagonal to 3D crystalline phase transition. Based on 



generic symmetry considerations of helical polymer ar- 
rays, we deduce that nascent crystalline groundstates 
fall largely into two categories: unfrustrated or weakly 
frustrated possessing hexagonal symmetry; and highly 
frustrated with rhombohedral symmetry. The minimal 
coupling allowed by symmetry between in-plane (bend- 
ing modes) and out-of-plane (torsion or sliding modes) 
deformation dictates the unusual elastic properties of 
DNA groundstates. In particular, for the case of highly 
frustrated groundstates, we find that the onset of order 
along the out-of-plane backbone direction dramatically 
alters the in-plane response of the array, leading to the 
emergence of a set of soft elastic modes in the crystal. 
Remarkably, the elastic response is highly resistant to 
those lattice distortions that transform crystalline pack- 
ing from one frustrated groundstate to another. 



I. INTERACTIONS AND GROUNDSTATES OF 
HELICES 

We begin by considering the interactions between 
aligned, single-helical molecules and demonstrate that 
crystalline groundstates of helical molecules on closed- 
packed lattices belong to two classes, unfrustrated and 
frustrated. In their most general form, the interactions 
between a pair of single helices, i and j, can be written 
as, 
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where U\/ are coefficients and 



and 



(1) 



are defined. 



Consider a horizontal plane intersecting the ith helix at 
some height, z = say [see Figure [1] (a)]. The azimuthal 
angle between vector connecting the central axis of the 
molecule to the spiral at that height and some reference 
direction (here y) defines fa. Note that due to helical 
symmetry of the molecule, fa can be changed by sim- 
ple rotation or translation of the molecule. This special 
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FIG. 1: Following ref. || the orientation of helical molecules 
is defined by projecting the helical groove location onto a bi- 
secting plane: (a) side view and (b) top view. The geometry of 
orientation for double-helical molecules is shown in (c), where 
4> g is the angular width of the minor groove. 



symmetry has been shown to the possibility of a peculiar, 
"screw liquid crystalline" phase of helical molecules § . 

In general, appearing in eq. ([T]) depends on de- 
tails of the inter- molecular forces present, helix separa- 
tion, and solution conditions such as concentrations of 
depletants or various ions species. Due to the finite- 
size molecule subunits of any helical molecules - such 
as, bases or G-actin monomers for the respective cases 
of DNA and F-actin - we can assume that series | uj? | 
is monotonically decreasing. This was shown explicitly 
for the case of the electrostatic interactions between uni- 
formly charged helices As such, we can assume that 
the (^-dependence of interactions are dominated by first 
few terms of eq. (flj. Below, we will focus on the terms 
n = 1 and n = 2, as these seem sufficient for describing 
observed DNA packings. 

Here, we avoid making specific assumptions about the 
nature of microscopic interactions. Progress can never- 
theless be made by considering two simple cases: repul- 
sive and "sticky" (or attractive) helices. Considering the 
geometry of Fig[T]it is clear that when <pi — (j>j = (2m+l)7r 
(to is any integer) the strands of neighboring helices are 
in close contact. For the case of repulsive (sticky) helices 
we therefore expect this close-contact orientation to cor- 
respond to a local maximum (minimum) of U s h- On these 
grounds we should expect that for the case of repulsive 
interactions, < when n is odd and U^' > for 

even n. Again, this expectation is borne out by electro- 
static calculations of like-charged helices. Of course, the 
opposite - ujfa > for rt odd and uj£' < for even n - 
would be true for sticky helices. 

For double-helical molecules like DNA, we can write 
the interaction in terms of inter-helix registry by sum- 
ming over the interactions between single strands. This 
yields an expression similar to eq. |T]), with single- 
helix coefficients replaced with double-helix coefficients, 
ujfa = 2J7 s k cos 2 (</> g /2), where 4> g is the azimuthal sep- 
aration between helical strands [see Figure [1] (c)]. As re- 



quired by symmetry, we see that for perfectly symmetric 
double-helical molecules, such as F-actin, with <j) g = n, 

the odd n terms of vanish. B-form DNA with a 
minor groove width of roughly of (j) g ~ 0.7tt — 0.87T, re- 
tains the n = 1 term, although with considerably dimin- 
ished amplitude. This, along with considerations of the 
foregoing paragraphs, suggests that we may capture the 
essential features of DNA packing with the simplified in- 
teraction potential, 

Udna ^ -Ui cos[0i - tj>j] + U 2 cos [2((j>i - #,•)] , (2) 

where U± and U2 are positive coefficients and we neglect 
the constant n = term. 

Focusing, in part, on the detailed dependence of Ui and 
U2 electrostatic forces present for DNA in aqueous solu- 
tion, Kornyshev and coworkers have demonstrated that 
interhelical potentials such as described by eq (HJ) lead to 
complex range of crystalline groundstates [TH EH, EH EH • 
For the case that U\ ^> U2, interactions are effec- 
tively ferromagnetic, XY-type and unfrustrated, leading 
to ground states with uniform azimuthal orientation. In 
the opposite limit, when U2 U\, interactions are frus- 
trated on close-packed, hexagonal lattices. Defining the 
angle = 2</>j, the second term in eq. @ is essentially 
an anti- ferromagnetic XY-type coupling in which 
cannot be mutually optimized be three pairs of adjacent 
molecules. Although this term prefers $i — $j ; = ir for 
all neighboring pairs, the best that can be accomplished 
on the hexagonal lattice with are configurations with <E>i 
changing by increments of ±27r/3 around the triangular 
plaquettes [14J. Hence, the frustrated groundstate breaks 
the full six-fold rotational symmetry of the hexagonal lat- 
tice, retaining only a three-fold rotation symmetry [see 
Figure [3 (a)]. Careful experimental studies DNA conden- 
sates [4, 5] reveal that close to the transition, crystalline 
structures maintain the hexagonal ordering of the lower 
density columnar-hexagonal phase. Only at higher densi- 
ties does the hexagonal order give way to an orthorhom- 
bic groundstate [Fig. (b)]. In order to resolve this 
structural transition in detail it is necessary to invoke a 
specific microscopic description of separation-dependence 
of interaction parameters, U\ and U2, as in the electro- 
static theory of refs. 0, EH, El El • 

An analysis of the groundstates of eq. © for helical 
molecules on the hexagonal lattice [ll[ reveals that the 
unfrustrated groundstate persists for Z7i > 4U2- Out- 
side of this range, the frustrated groundstates consist 
of three staggered sublattices-labeled as xa , xb and xc 
respectively-as shown in Figure [21 If helices on sublat- 
tice A have orientation, </>±(xa) = </>o, then those on the 
respective B and C sublattices are given by </>±(xb) = 
4>o ± S(p and </>±(xc) = 4>o T f>4>, where the value of 5<p 
varies continuously from the unfrustrated state, S(f> = 0, 
when Ui = AU 2 to the fully-frustrated state, 6(f) = ir/3 
when U2 S> Ui. Here, ± refer to the two degenerate 
groundstates, which differ by a rotation in the plane by 
2tt/6. 
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II. CRYSTALLIZATION OF UNFRUSTRATED 
GROUNDSTATES 

Having considered the groundstate structure of helical 
arrays, we now turn to construct a proper description of 
transition from the 2D columnar-hexagonal phase to a 
higher density, 3D frustrated or unfrustrated crystalline 
phase. We begin by constructing an order-parameter 
theory for the transition to the simpler unfrustrated 
groundstate, where helices have a uniform orientation 
(i.e. <f>i — const, for all i). 

In the absence of the out-of-plane order, the columnar 
hexagonal phase is described by the displacement field 
u(r) that describes fluctuations of the polymer backbone 
in the xy plane of hexagonal order. Here, r = x + z, 
where x labels the in-plane position of the chain in the 
hexagonal lattice and z describes the position along the 
polymer backbone. The (coarse-grained) elastic energy 
describing fluctuations of this two-dimensionally ordered 
phase can be written as 



(b) 



Hx = 



2pu% + Xu 2 kk + K 3 (d 2 z uf 



(3) 



FIG. 2: (a) shows the symmetry of the rhombohedral ground- 
state of helical packing. The sublattices A, B and C refer 
relative rotational displacements of helices. In the case of 
extreme frustration, sublattices B and C are respectively dis- 
placed ±7r/3 (or =p7r/3) relative to A, where p is the helical 
pitch, (b) shows are alternative orthorhombic groundstate 
which alters the underlying hexagonal lattice of lower density 
columnar phase. The relative displacement of sublattices A 
and B are dictated by details of the interhelical potential. 



Clearly, groundstate properties helical molecular as- 
semblies are most sensitive to the precise values of U\ 
and U2 for DNA and other biopolymers, and these co- 
efficients themselves will have a complex dependence on 
microscopic parameters resulting from a combination of 
intermolecular forces. Even without a full knowledge of 
the inter-molecule forces at play in biopolymer systems, 
we can make progress by noting that groundstates of heli- 
cal molecules will larg ely fall into two classes, frustrated 
and unfrustrated [2lfl . The symmetry of the former is 
hexagonal, while the latter possess only a rhombohedral 
symmetry. In the remainder of the article we focus on the 
emergent elastic properties of these groundstates based 
on the broken symmetry of crystalline groundstates. A 
similar symmetry-based approach was used by Lorman, 
Podgornik and Zeks to explore the thermodynamics of 
the transitions between low-density hexatic phases and 
high-density orthorhombic phases of DNA [15|]. Unlike 
that scenario, we find below that it is not, in general, 
possible accomplish that transition continuously if the 
system first passes through an intermediate rhombohe- 
dral packing. 



Here, i, j and k refer only to in-plane directions and 
Uij = (diUj +djUi)/2 is the symmetric, 2D elastic strain. 
The first and second terms describe distortions of the 
hexagonal packing and the third term describes the en- 
ergy cost for fluctuations that bend the polymers. In 
principle, this energy may also include additional penal- 
ties for so-called twist and splay deformations [171 ]. 

In the unfrustrated groundstate, the onset of order 
along the polymer backbone direction can be described 
by the development of a non-zero amplitude of a smectic- 
like order parameter, ^(r) p7| . Specifically, we can define 
the order parameter, which represents modulations of 
the material density, p(r) along the z direction, through 
p(r) ~ Re[l+V'( r ) elG2: ]pHex(r), where pHex(r) is the den- 
sity field of the 2D array of polymers. The special sym- 
metry of helical molecules with pitch, p, is captured by 
the following transformation properties of the order pa- 



rameter. Under translation along z by Az, ip — > e 



;GAz 



1>, 



where G — 2tt /p. And under rotation around I by A#, 



AA8 



ip. This ensures 1) that the order parameter 



has the appropriate periodicity and 2) that a transla- 
tion followed by the appropriate rotation leaves a helical 
groundstate unchanged. 

Since the order parameter of the unfrustrated ground- 
state is essentially that of a 3D XY model, the following 
Landau free energy describes the onset of a non-zero am- 
plitude of tp( r ) a t the columnar hexagonal to crystalline 
transition, 



H 



I/- 



(4) 



In addition to H|| ps j| we must also consider gradient terms 
that are consistent with the symmetry of both phases. 
On it's own eq. ^ has the symmetry properties of 
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a 2D hexagonal array of chains, namely, the energy is 
invariant under rotations around z by 2ir/6. The un- 
frustrated groundstate maintains this in-plane symmetry, 
only breaking translational symmetry along the back- 
bone direction. To second order in derivatives we have 
the following contributions to the free energy of the array, 

+bu kk Re[iipd z il)*} + \ip\ 2 (28p u kk + SX u 2 k + 2Sp w? )}. 

Here, cm, c±, b, Sp, Sp and SX are coefficients. Note 
the appearance of the gauge coupling between in-plane 
gradients of ip and z derivatives of backbone displace- 
ment in the second term of (O- This is identical to the 
rotationally-invariant coupling between gradients of the 
order parameter and the director for smectic phases [TtJ , 
as d z u is the in-plane component of the backbone tan- 
gent vector. It is this coupling which is responsible for 
the resistance of the crystal to_ shear deformations that 
slide molecules past another (l8| . 

The specific critical properties of the phase transition 
from 2D columnar tho 3D solid order are unusually com- 
plex, owing to the gauge coupling between in-plane dis- 
placement modes and fluctuations of the smectic-like or- 
der parameter. For the purposes of the present study 
we focus on emergent properties outside of the critical 
region. In the crystalline state (above a critical density) 
if) takes on a non-zero amplitude, \ipo\ = (— r/it') 1 / 2 , for 
r < 0, where u' = u — 5p/A{X + p). To examine fluctu- 
ations of the groundstate, we consider phase variations 
of the order parameter, ip(r) — \il)o\e l ^ r \ which repre- 
sent torsional or translational displacements of the helical 
macromolecules. Defining C\\ = C|||V»o| 2 , C± = c±\ip \ 2 , 
2B = 6|^ | 2 , A' = A + \ip \ 2 5X and p! = p + \^ \ 2 Sp we 
have the following effective Hamiltonian for in-plane (u) 
and out-of-plane (0) fluctuations, 



n ef f[u,^}=-J d'x^id^Y + C ± (V^ + Gd z u) 2 (6) 
+2Bu kk d z <P + 2p' u% + X'u 2 kk + K 3 (d 2 u) 2 

Here, we will focus on the renormalized elastic behavior 
of the in-planc displacement modes due to the coupling 
between the out-of-plane modes appearing in the crys- 
talline phase. Integrating out the smooth (Gaussian) 
fluctuations of the 4> field in it is straightforward 
to derive an effective Hamiltonian (in Fourier space) for 
the u fluctuations in the crystalline phase, 

K ff \u] = l I S^[Sz\^)\u L (k)\ 2 +S^(k)\u T (k)\ 2 }, 



and 



^ff^-2j W) 



where = kj^(k^ • u)/^^, uy = u u^, 

S- L \K)={X' + 2p')kl+K 3 k 4 z 

C\\C ± G 2 k A z - B(2C\G + B)k 2 z k\ 
+ Cnk 2 + C^k 2 



(7) 



(8) 



S^(k) =p'k] 



Kokt 



C±G z k 



2;„2 



(9) 



The transverse modes, ur(k), are generically stable and 
behave as simple phonons in an anisotropic 3D crystal, 
while the fluctuation spectrum of the longitudinal modes 
is considerably more anisotropic owing to the coupling to 
out-of-plane modes. 



III. CRYSTALLIZATION OF FRUSTRATED 
GROUNDSTATES 

Having analyzed the simpler unfrustrated case, we now 
construct the Landau-Gizburg approach to the columnar 
to crystalline transition when the groundstate is highly 
frustrated. Here, the groundstates have the staggered 
structure shown in Figure [2] (a), with helical molecules 
occupying 3 sublattices. The symmetry of the ground- 
states is encoded by the phases </>±(x) discussed above, 
so that the order parameter for the frustrated states can 
now be written as, 



V>(r) = V+(r)e^+ (x) + V-(r)e^- (x) 



(10) 



where again ± corresponds to crystalline ordering in ci- 
ther of the two degenerate possibilites. Hence, for highly 
frustrated helical groundstates we have a 2-component 
complex vector order parameter (V>+, ip-)- 

As both groundstates have rhombohcdral symmetry, 
they are invariant under rotations in the plane by 2ir/3, 
but they are not six-fold symmetric. Consider a rotation 
in the plane of 2tt/6 around an axis passing through the 
A sublattice. From Fig. ^ it is clear this operation takes 
A sites to A sites and interchanges B and C sites. Since 
0±( x b) = — <^±(xc), under roations by 2ir/6 in the plane, 
(■0+, V 1 -) ~~ * (V'-) ( U P to an overall phase shift due to 
helicity of the molecules). In constructing the Landau- 
Ginzburg free energy for the transition, all terms must 
respect the six-fold rotational symmetry of the higher- 
symmetry, columnar hexagonal phase. To fourth order 
in 1^1 4 the most general Landau free energy describing 
this type of U(l) x Z2 symmetry breaking is jTgJ, 



n 



1 



d 3 [r(\^ 



+ 2(l^+l 2 + IV;-l 2 ) 2 + ^ + l 2 IV'-l 2 ]. (11) 



Note that r H\ x p ± \ is invariant under interchange of + and 
— components (rotations by 27r/6). Spontaneous sym- 
metry breaking is described by the case r < 0, u > and 
v > —2u [v = —2u,r = is a tricritical point). In this 
case, the order parameter assumes a non-zero amplitude 
in either the + or — component, while the other compo- 
nent goes to zero. Here, this broken discrete symmetry 
indicates the broken six-fold rotation symmetry of the 
hexagonal columnar phase. 
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For the case of the frustrated groundstate, the same 
gradient terms appearing in eq. ([5]) for the unfrustrated 
case will also appear in TL^j for the frustrated case, except 
that each term with tp in the former case is replaced 
by term with both tp + and maintaining + and — 
exchange symmetry. In addition to the gradient terms 
consistent with hexagonal symmetry, we are obliged to 
consider coupling terms that are consistent with the lower 
three-fold symmetry frustrated groundstates. This term 
is constructed from products of the vectors, 



,± _ 



Re[#±(Vx+iG5 z u)^y, 



(12) 



and Uij which transform as a vector and tensor, respec- 
tively, in the xy plane. Following ref. 20], it can be 
shown that the three-fold rotational symmetry of the 
rhombahedral state allows for two possible free-energy 
contributions, 



"mi 



]+A±[(t 



-Uyy)v±+2U XV vt] , 

(13) 

where X x and are coefficients. One of these two terms 
must be eliminated due to the symmetry properties of 
frustrated groundstates of helical molecules, which are 
invarient under the combined operation y — > — y and 
z — > — z (see the in-plane mirror symmetry in Fig. [2]). 
Note, that the second reflection is required by the chi- 
ral symmetry of the molecules. Since vf~ — —vf' under 
z — > — z, this requires that A^ = as this term changes 
sign under the mirror reflections, leaving only A^ ^ 0. 
Of course, for a crystal of achiral molecules, the opposite 
scenario would be true (A^ ^ and A^ =0). 

To be clear, we have constructed elastic terms, 
^ \iuxx — Uy y )v y ^ + 2u xy v x ] that respect the three-fold 
symmetry of the frustrated crystalline states, tp±. Under 
rotations by 2w/6, this term changes sign. Since the + 
and — states are exchanged under this operation, choos- 
ing A + = — A~ = Xf yields a six-fold rotationally invari- 
ant contribution to the free energy, 



5H\ 



A, 



d A x 



( V V - V V ) ( U XX -Uyy)+2(V+-V x )l 



(14) 

As before we examine the effects of fluctuations of the 
out-of-plane and in-plane degrees of freedom in the bro- 
ken symmetry crystalline state, by considering phase 
variations of the order parameter. If we take as our 
broken symmetry state, say, - 0+( r ) = iV'ole 1 ^ 1 ' and 
■0-(r) = 0, we then have the following additional elastic 
contribution to the fluctuation Hamiltonian, TL e ff[<f), u], 
for frustrated groundstates, 



5H e ff[(f>,u}=Af / d 3 x (d v 4> + Gd z u y )(u xx - u yy ) 



+2(d x (t) + Gd z u x )u xl 



(15) 



where A/ = A/|?/>o| 2 - 

The appearance of this new coupling between (j> and 
u modes in eq. dT51) reflects a highly anisotropic elastic 




FIG. 3: The wavevectors and polarizations of soft elastic 
modes of the frustration rhombohedral state are shown in (a) 
and (b) respectively. In (c), the geometry of the soft simple 
shear modes is shown on the left, while the elongation shear 
that carries the rhombohedral state [Fig. [2] (a)] to the higher 
density orthorhombic state [Fig. [2] (b)] . 



response of the frustrated crystalline state. Integrating 
out the out-of-plane {4>) modes yields expressions for the 
in-plane response of the system even more complex than 
eqs. (HI)-®- We can nevertheless get a sense of the elastic 
response of the frustrated states by focusing on a simple 
limits of the full elastic response as defined by, 



i r rl 3 k « , 



(16) 



where i and j refer to independent directions in the plane. 
The components of (k) can be extracted from eqs. ([6]) 
and (fT5|) . Specifically, we focus on the k z = terms in the 
limit of incompressibility, ■ u(k) = 0. For this special 
case, (k) is diagonal in the basis, Mi(k) and «2(k) 
which can be written in terms of x and y displacements 
as, 



/, % 2k x k y (k x k ) 
Ul ( k ) = 1,2 M k ) H ,.2 u v\*)> 



«a(k) = 



kl 
(k\ — ky) 



u x (k) 



2/C x ky 



u„(k), 



with eigenvalues, 



5 1 - 1 1 (kx,fe = 0)= ^'-^Ijfel, 
S 22 1 (kx,fe z = 0) = M 'fcx' 



(17) 
(18) 



(19) 
(20) 



Here, we find an especially soft compliance along the 
ui(k) directions. The shift of (k.±,k z = 0) due to 
the coupling to out-of-plane 4> modes is proportional to 
t/>o| 2 j indicating the splitting between the response to the 
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Ui and U2 modes grows with increasing crystallinity. The 
incompressibility constraint, kj^ - u = implies a relation 
between the principle directions of in-plane deformation. 
This constraint only allows for deformations with Ui ^ 
while U2 = when k y = or when 3fc 2 = ky. Hence, 
the wavevectors of these soft modes lies along the crystal 
axes of the underlying hexagonal lattice (see Figure [3]). 
From eq. (TP?)) we find that the polarization of the soft 
modes is perpendicular to its direction in k-space. That 
is, upon crystallization to frustrated groundstate, the he- 
lical array is most easily deformed by simple shear with 
the shear direction perpendicular the hexagonal axes. 

It must be emphasized that special compliance to sim- 
ple shearing of the underlying hexagonal lattice of the 
frustrated state is a necessary consequence of the discrete 
breaking of 6-fold symmetry and does not rely on details 
of the microscopic interactions between DNA molecules. 
This analysis relies on only 1) the 3-fold rotational sym- 
metry of the frustrated crystalline groundstate 2) the 
underlying 6-fold symmetry of columnar phases and 3) 
the lack of mirror symmetry of the constituent chiral 
molecules 22]. A direct signature of this anisotropic 
phonon dispersion appears in the diffuse thermal scatter- 
ing surrounding Bragg peaks in the k z = plane. Near 
to a Bragg spot at reciprocal vector, G, this scattering 
is encoded in (|G • u(Ak)| 2 ), where Ak is the wavevec- 
tor deviation from G. This analysis predicts pronounced 
thermal scattering (for Ak _!_ G) around Bragg peaks 
whose G lies along one of the six soft-mode polarizations 
(Fig. [3]), and reduced thermal scattering for G off those 
axes. 

What is perhaps most surprising is the relationship 
between the anisotropic elasticity of the DNA array and 
the deformations which carry the systems from one frus- 
trated groundstate to another. It is clear from Fig. [2]that 
the rhombohedral groundstate is most directly trans- 
formed into the higher density orthorhombic groundstate 
by an elongational shear, i.e. u y = ~ey and u x = ex [see 
Fig. [3] (c)]. According the linear theory, such modes 
lead to particularly large elastic energy costs. The ratio 
of the elastic cost of this mode relative to simple shear 
(i.e. u y — ex and u x — 0) is 9(2 — a)/(l — a) > 18, 
where a = A|/(//Cj_). While the linear-elastic theory of 
the frustrated rhombohedral DNA crystal does not pro- 
vide a full description of the structural transition to or- 



thorhombic groundstate, favored at higher densities [|[, it 
nonetheless makes certain generic predictions regarding 
the dynamics of this phase transition. Specifically, we ex- 
pect the rhombohedral DNA crystal must pass through 
some intermediate state associated with a simple-shear 
deformation before reaching the orthorhombic state since 
the crystal is mechanically least stable to these modes. 
Interestingly, the orthorhombic state is not accessible 
by pure simple shear alone, and secondary distortion is 
required to accomplish this structural transition. This 
should be distinguished from scenario where the system 
bypasses the rhombohedral state and may pass from hex- 
atic to orthorhombic crystal through series of continu- 
ous transitions [X5j | . In light of this prediction, it would 
interesting to resolve experimentally the structural in- 
termediates which accommodate the observed change in 
symmetry between these states of DNA packing. 



IV. SUMMARY 

Relying on symmetry considerations alone, we have 
constructed an elastic theory of DNA crystals, which are 
inherently frustrated by the underlying helical structure 
of the constituent molecules. As the system passes from a 
hexagonally-ordered liquid crystalline state to a 3D crys- 
tal at higher density, the interactions between double- 
helical molecules are, in general, quite frustrated leading 
to spontaneously broken discrete symmetry associated 
with the lower 3-fold symmetry of the groundstate. A 
general consequence of the of this broken symmetry is 
a highly anisotropic response to deformations of in-plane 
order, with a pronounced resistance to elongational shear 
modes associated with structural transitions to higher 
density groundstates. 
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